Spin excitations in fluctuating stripe phases of doped cuprate superconductors 
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Using a phenomenological lattice model of coupled spin and charge modes, we determine the spin 
susceptibility in the presence of fluctuating stripe charge order. We assume the charge fluctuations 
to be slow compared to those of the spins, and combine Monte Carlo simulations for the charge order 
parameter with exact diagonalization of the spin sector. Our calculations unify the spin dynamics 
of both static and fluctuating stripe phases and support the notion of a universal spin excitation 
spectrum in doped cuprate superconductors. 



A key challenge in the field of high-T c superconduc- 
tivity is to separate universal from non-universal proper- 
ties. For spin fluctuations, believed to be the glue that 
binds the Cooper pairs, this issue is controversial: early 
neutron scattering experiments had established the ex- 
istence of a "resonance peak" , corresponding to a spin 
collective mode at the antiferromagnetic wavevector, for 
certain cuprate families [1-3], while in others stripe- 
like spin and charge modulations were detected [4-7]. 
(Signatures of charge order, likely pinned by impuri- 
ties, have been observed also in scanning tunneling mi- 
croscopy (STM) experiments on Z^S^CaC^Os+a [8] 
and Ca2- a; Na a ;Cu02Cl2 [9].) Recent experiments have 
mapped out the spin excitations in various cuprates over 
a large range of energies [10-14], with remarkable results: 
(i) In stripe-ordered La 1 5/ 8 Ba 1 / 8 Cu04 (LBCO) two ex- 
citation branches have been found, with the high-energy 
branch above a "resonance" well described by the spec- 
trum of a gapped spin ladder [10]. (ii) YBa2Cu306+a 
shows incommensurate excitations below the resonance 
energy [11-13]. These results point toward a low- 
temperature spin excitation spectrum being universal 
among the cuprate families at intermediate energies [15], 
namely an "hour-glass" spectrum with a high-intensity 
peak at wavevector (w, tt) and both downward and up- 
ward dispersing branches of excitations. 

However, a unified theoretical description for the spin 
dynamics is lacking. The spin excitations in LBCO [10] 
appear to be well explained within a model of static 
stripes [16-18], where weak magnetic order exists on top 
of a bond-ordered (i.e. dimerized) background. In con- 
trast, neutron scattering in YBa2Cu306+a was modelled 
using RPA-type calculations [19, 20] which, however, rely 
on details of the band structure and are not able to de- 
scribe ordered states. A key question is whether the neu- 
tron scattering data on YBa2Cu3 06+5 can be understood 
in a stripe picture as well - as no static order has been 
detected, stripes have to be fluctuating in space and time 
here. Controversial experimental viewpoints on this have 
been put forward [6, 11, 15]. 

The purpose of this letter is to show that fluctuat- 
ing stripes [21-23] lead to spin excitations very similar 



to those observed in the experiments, thus providing a 
unified account of the collective mode dynamics in the 
cuprates. While large, weakly fluctuating, stripe do- 
mains are consistent with the results in L^-zSr-rCuG^ 
and La2- a ;Ba a ;Cu04, a mixture of stripe and checker- 
board structures is required to explain the experimental 

data on YBa 2 Cu3 6.85 [H]- 

Lattice order parameter theory. We employ a phe- 
nomenological model of coupled spin and charge fluc- 
tuations [24] where the spin incommensurabilities are 
driven by inhomogenieties in the charge sector; this 
is supported, e.g., by experiments on stripe-ordered 
La2_ a; Sr a ;Cu04, where the charge order sets in at a higher 
temperature than the spin order. On a microscopic scale, 
the influence of the charge order on the spin sector can be 
understood as spatial modulations of both spin densities 
and magnetic couplings [16, 17]. (Additional collective 
degrees of freedom with zero wavevector, e.g., pairing 
fluctuations, will not qualitatively modify our results.) 
The goal of our work is to describe well-defined collective 
modes, hence we neglect the continuum of single-particle 
excitations and the associated collective mode damping. 

The action of our Landau theory has the form S = 
S v + S^, +S v ^p, where S V (S^) describe the spin (charge) 
fluctuations, and S v ^p couples the two. We assume a 
dominant antiferromagnetic interaction, and so employ 
a lattice tp 4 theory for the spin fluctuations at the com- 
mensurate wavevector Q — (7r,7r), 

S <P = J dT ^ t( 9 ^') 2 + s $] + XX^' -$3'f + S i 

3 (33') 

with 1S4 being the quartic self-interaction term. The real 
order parameter (p 3 - and the spins Sj on the sites j of the 
square lattice are related through Sj oc e l< ^' rj 0j. 

Turning to the charge sector, we note that micro- 
scopic calculations have indicated a tendency towards 
states with stripe-like charge ordering [25, 26], but states 
with two-dimensional (2d) "checkerboard" modulations 
closely compete in energy [27-29]. We employ two com- 
plex order parameter fields ip x ^ y (r,T) which measure 
the amplitude of horizontal and vertical stripe order at 



wavevectors K x ^ y . Checkerboard order then implies both 
ip x and ip y non-zero. In a situation with fluctuating 
charge order, the balance between stripes and checker- 
board (which depends on microscopic details [27] ) is con- 
trolled by a repulsion or attraction between tp x and ip y . 
The complex phase of i/) Xjy represents the sliding degree 
of freedom of the density wave and distinguishes between 
bond- and site-centered stripes. In our simulations we 
concentrate on K x — (n/2,0) and K y = (0, 7r/2), i.e., a 
charge modulation period of 4 lattice spacings; modula- 
tions at these wavevectors have been observed both in 
neutron scattering [4-7] and STM [8, 9], in particular 
near doping 1/8 where stripe order is most robust. The 
real field Q x {r) = Reip x (r)e iK,c ' r (similarly for Q y ) mea- 
sures the modulation of both the charge density (for r 
on sites) and bond order (i.e., kinetic energy or pairing 
amplitude, for r on bonds). We choose signs such that 
5p{rj) = Q x + Q y is the deviation of the local hole density 
from its spatial average. The all-important couplings be- 
tween spin and charge fluctuations have to be of the form 
XQip 2 due to the underlying SU(2) symmetry. Guided by 
the lattice models [16, 17] we choose [24] 
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Ai > implements the correlation between the on-site 
charge density and the amplitude of the spin fluctuations, 
while A2-4 ensure that the effective first- and second- 
neighbor exchange constants modulate along with the 
bond order; the antiphase domain wall properties of the 
stripes [7, 25] are reflected in the positive sign of A 2 ,3. 

For constant t/j XiV the action S v + is a theory for 
magnetic modes in a background of static charge order. 
For sufficiently large A couplings, the minimum energy of 
the tp fluctuations will be shifted away from (jr, it) to the 
incommensurate wavevector dictated by the charge order 
(this is a non-perturbative effect!), with the spin order 
remaining collinear. Results [24] for the spin fluctuation 
spectrum in the presence of static stripes are in excellent 
agreement with results on LBCO [10]. 

Fluctuating charge order. We now turn to <S^: For 
slowly fluctuating charge order it is useful to think about 
snapshots of the charge configuration. Defining an 0(4) 
field ip = (ip x , ip y ) we can discuss physically distinct spa- 
tial fluctuations: (i) Fluctuations of the complex phases 
of ip XtV are stripe or checkerboard dislocations, (ii) Fluc- 
tuations between areas of dominant tp x or tp y represent 
domain walls between horizontal and vertical stripes, (iii) 
Variations in |^| are amplitude fluctuations in the local 
charge order. Depending on the particular form of S^, 
these fluctuations will have different importance. Re- 
garding amplitude fluctuations two extreme cases come 





FIG. 1: (color online) Snapshots of the charge order param- 
eters \ip x ,y\ and the resulting charge modulation (Q x + Q y ), 
obtained from MC simulations for bond-centered fluctuating 
stripes on 32 2 sites, using a ip 6 action S^, see text, a) Strong 
repulsion between ip x and ip y , sharp domain walls, b) Weak 
repulsion, smooth domain walls with checkerboard structure, 
c) Weak attraction, fluctuating checkerboard order. 

to mind: (a) a standard ip 4 theory with a "soft" or- 
der parameter, which has rather large amplitude fluc- 
tuations, and (b) a "hard" order parameter theory with 
a fixed- length constraint, \ijj\ 2 — const. Although mi- 
croscopically amplitude fluctuations are present, existing 
approximate results for Hubbard or t-J models are in- 
conclusive with regard to their importance. Experimen- 
tally, STM results [9] on Ca 2 -xNa :c Cu02Cl2 indicate a 
spatially disordered arrangement of stripe segments and 
more 2d "tiles" , with the amplitude of these local mod- 
ulations fluctuating rather little. 

In our simulations, we have employed various forms 
for the charge action <S^, with different amounts of am- 
plitude fluctuations. Most useful is a ?/> 4 -typc theory for 
the 0(4) field ip, supplemented by a positive ip 6 term: 
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with ip 2 = \ip x \ 2 +\ip y \ 2 . A combination of u\ < and u 2 > 
suppresses amplitude fluctuations of tp. For c\ x = c\ y , 
c-2x = C2y, s x = s y , andi> = w = 0, the action has 0(4) sym- 
metry. The w term selects between bond-centered and 
site-centered stripes. The important quartic ^[V'zPlV'i/l 2 
term regulates the repulsion or attraction between hori- 
zontal and vertical stripes, i.e., it determines whether the 
character of the order will be one-dimensional (stripe, for 
v > 0) or two-dimensional (checkerboard, for v < 0). 
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FIG. 2: (color online) Dynamic susceptibility \'{q,oj) 
for bond-centered fluctuating stripes [31] on 40 2 sites. 
Left/Middle: cuts at a constant energy, slightly below/above 
the resonance energy, E las , as function of momentum. Right: 
cuts along (q x , 7r) as function of q x and energy, showing the 
universal "hour-glass" spectrum. The couplings are Ai = 
A3/2 = S.Eres/l'i/'ltyp, A2 = A4 = 0. a) Strong repulsion be- 
tween tp x and ipy, correlation length £ « 30. b) Weak repul- 
sion, £ « 20. c) Weak attraction, £ ~ 20. 

To simplify the treatment of S v +S^, +S V ip, we assume 
that fluctuations in the charge sector are slow compared 
to those in the spin sector, and we neglect the feedback 
of the spins on the charges. This leads to an adiabatic 
(Born-Oppcnhcimer) approximation for the coupled dy- 
namics, and allows to treat the charge fluctuations by 
classical lattice Monte Carlo (MC) simulations. For each 
configuration of the ip x ,y, the remaining theory S v 
(at the Gaussian level, £4 = 0) is quadratic in the <p 
fields and can be diagonalized on lattices up to 64 2 sites. 
(Neglecting 1S4 is justified in spin-disordered phases.) We 
employ a standard Metropolis algorithm with single-site 
updates at a finite effective temperature (T = 1) to sim- 
ulate <S^ in a regime where the correlation length £ is 
between 5 and 50 lattice spacings [30]. The spin suscep- 
tibility %"(g, w) is obtained by averaging its value over 
typically 20 MC charge configurations, with 10 5 - 10 6 
MC steps between two measurements. 

Numerical results. Typical snapshots of the two charge 
order parameters tp x ,y an d the resulting charge configu- 
ration, for different values of the stripe interaction v, are 
shown in Fig. 1. Let us now discuss our results for the dy- 
namic spin susceptibility, \" (q, lo), as measured in inelas- 
tic neutron scattering. Starting from ordered stripes [24] , 
we found that amplitude fluctuations of tp rather quickly 
destroy the incommensurate spin response; for a stan- 
dard ip 4 theory this happens already within the ordered 



phase. In contrast, the spin sector turns out to be less 
sensitive to phase fluctuations of the stripe order. Inter- 
preted microscopically, this means that incommensurate 
spin response requires well-formed stripe segments with 
a length of order 10 lattice spacings. We have therefore 
focussed on versions of with small amplitude fluctua- 
tions, and carried out large-scale simulations for various 
couplings and correlation lengths of the il> x ,y 

Results for the dynamic susceptibility, corresponding 
to the situations in Fig. 1, are shown in Fig. 2. The 
right panels show that a common feature of all spec- 
tra is an "hour-glass" (or "X-shaped" ) spectrum with a 
strong "resonance" peak at (tt, tt) and energy E les . The 
downward dispersing lower branch is most pronounced 
for well-defined stripes, i.e., large £; it is progressively 
smeared out with decreasing £ (Fig. 3). Very recent 
neutron scattering results [32] indicate a "Y-shaped" re- 
sponse in the pseudogap state above T c in underdoped 
YBa2Cu306.6, with little dispersion at low energies - 
these data seem to be consistent with fluctuating short- 
range stripe segments, Fig. 3c. 

We have studied in detail the crossover from a strictly 
"stripy" situation with a strong repulsion between ip x and 
ip y (Figs, la, 2a, 3) to a 2d checkerboard regime with at- 
traction between tp x and ip y (Figs, lc, 2c). Significant dif- 
ferences occur in the lower branch (left panels of Fig. 2). 
For large domains of horizontal or vertical stripe order, 
Fig. 2a, well-defined peaks occur along (q x ,n) and (tt, q x ), 
as observed in LBCO. Interestingly, with increasing vol- 
ume fraction of checkerboard domain walls, Fig. 2b, the 
neutron response is both smeared and enhanced along 
the g-space diagonals, resulting in a quasi-2d dispersion 
of the downward branch - this is strikingly similar to 
data on YBa2Cu306.85 [12]. Finally, in a checkerboard 
regime the low-energy spin excitations occur along the 
diagonals (Fig. 2c). Based on these results and STM 
data [9] we predict that such a neutron response should 
be observable in Ca2- x Na. I ;Cu02Cl2. 

The g-space structure of the upward dispersing branch 
(middle panels of Fig. 2) changes less from the stripe 
to the checkerboard regime; it is strongly anisotropic 
only in the stripy situation of Fig. 2a where it resembles 
the spectrum of two-leg ladders. Focussing on the right 
panels of Fig. 2, we further observe that larger checker- 




FIG. 3: (color online) Evolution of x"(<fj tJ ) with decreasing 
stripe correlation length £, for a situation with strong repul- 
sion between ip x and rp y , as in Fig. 2a. 



4 



(0,2rc) (7x,2tz) 




(0,0) (2it,0) 



FIG. 4: (color online) Dynamic susceptibility slightly below 
the resonance energy, for fluctuating stripes in the presence 
of an in-plane anisotropy. From a) to c) the anisotropy is 
increasing: the ratio of the gradient coefficients for tpx,y in 
the charge action is a) f .005, b) 1.01, c) 1.02. 

board regions tend to suppress the upper branch right 
above the resonance peak, it re-appears only at some- 
what higher energies. This is in remarkable agreement 
with neutron data on YBa2Cu30g.85 [14], and can be 
easily understood: For perfect checkerboard order the 
low and high-energy response are separated by a large 
gap [24] , and our simulations interpolate between stripes 
and checkerboard. (An alternative interpretation of the 
data of Ref. 14 within RPA is in Ref. [20].) 

To model detwinned YBa 2 Cu3 06.85 [11], it is neces- 
sary to include an in-plane anisotropy to account for the 
orthorhombic distortions. Assuming the anisotropy to 
be small, it will mainly influence the low-energy charge 
fluctuations. We have therefore carried out simulations 
where the ip x and ip y order parameters had different mass 
and/or different gradient terms (i.e. velocities). Sample 
results for the downward dispersing branch are shown in 
Fig. 4, which are in reasonable agreement with the data 
of Hinkov et al. [11]. The anisotropy decreases at higher 
energies (not shown). 

Let us note that our adiabatic approximation to treat 
Stp + Sip + cannot distinguish between slowly fluc- 
tuating and time- independent disordered stripes (e.g., a 
"stripe glass" pinned by impurities). Physically, these 
two situations will indeed yield very similar spin excita- 
tions; a clear-cut distinction will require a direct probe 
of the charge modes. 

Conclusions. We have determined the spin excitations 
in the presence of fluctuating stripe charge order. We 
obtain an incommensurate spin response provided that 
(i) charge order fluctuates predominantly in phase rather 
than in amplitude, and (ii) the charge correlation length 
is at least 10 lattice spacings. (Assuming a collective 
mode velocity of 50 meV this roughly translates into THz 
fluctuation frequencies.) In addition, we found that an 
increasing volume fraction of stripe domain walls with 
checkerboard structure leads to quasi-2d spin fluctuations 
as observed in YBa2Cu3 06.85- 

Our calculations thus support the notion of a universal 
spin excitation spectrum at intermediate energies in the 
cuprates, arising from stripe-like charge-density fluctua- 
tions. This brings us closer to a unified description of the 
collective excitations in the high-T c materials. 
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